Theory of electric field induced one-magnon resonance in cycloidal spin magnets 
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We propose a new mechanism to induce a novel one-magnon excitation by the electric component 
of light in cycloidal spin states, i. e. so called electromagnon process. We calculated optical spectra 
in the cycloidal spin structures as observed in multiferroic perovskite manganites -RMnOa where 
novel magnetic excitations induced by oscillating electric fields are observed. When symmetric spin- 
dependent electric polarizations are introduced, we have light absorptions at terahertz frequencies 
with one- and two-magnon excitations driven by the electric component of light. Our results show 
that some parts of optical spectra observed experimentally at terahertz frequencies are one-magnon 
excitation absorptions. 

PACS numbers: 75.80.-|-q, 75.40.Gb, 75.30.Ds, 76.50.-|-g 



MultifFeroics, where fcrromagnctism and fcrroelectric- 
ity coexist, have attracted both experimental and the- 
oretical interests d, 0]. Recently, ferroelectric per- 
ovskite manganites i^MnOa {R = Tb, Dy, Gd, and oth- 
ers) are especially paid attention to due to colossal mag- 
netoelectric effects Q , since ferroelectricities are strongly 
associated with their cycloidal spin orders, and as a re- 
sult, they have strong magnctoclcctric couplings. The- 
oretically, such a strong coupling between the ferroelec- 
tricity and the cycloidal spin state can be understood by 
considering a spin-current induced polarization described 
by Pa oc x {Si x Sj) @, Once the cycloidal spin 

state is realized, an electric polarization breaks out whose 
direction is perpendicular to the propagation and the he- 
licity vectors of the cycloidal spin structure. One of the 
experimental evidences to confirm the spin-current sce- 
nario is that, in /JMnOa, the direction of polarizations 
can be controlled by changing the direction of the cy- 
cloidal plane with external magnetic fields Q]- For ex- 
ample, in TbMnOs the polarization flops from P\\c to 
P\\a, by a magnetic field induced spin flop from he to ab 
cycloidal states. 

One of the hot topics in these multiffcroics materials 
/JMnOa is electric field induced spin excitations, or elec- 
tromagnon, observed _in _aii optical spectroscopy at tera- 
hertz frequencies 
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neering work by Pimenov et at, they observe magnetic 
absorptions in TbMnOs and GdMnOs for i;||a Espe- 
cially, absorption around 2 meV in TbMnOs is expected 
to be strongly related to ferroelectricity induced in the 
spin-current model. In fact, the spin-current model pre- 
dicts that the electric component of light perpendicular 
to the cycloidal spin plane, i.e. E\\a in TbMnOa, can ex- 
cite a one-magnon excitation [l5| . Thus, excitations ob- 
served in TbMnOs were considered as such one-magnon 
excitations. 

However, further detailed investigations revealed that 



terahertz absorptions in /JMnOa could not simply be un- 
derstood by the theories proposed so far. Main unclear 
feature is concern with a selection rule in the spiral or- 
dered phase, which does not depend on R site ions, tem- 
perature, external magnetic fields, and hence directions 
of spins. This is in contradiction with the prediction of 
the spin-current model E Jl P. For example, at the low- 
est temperature region of the spiral spin state in zero 
magnetic field, P\\c in TbMnOa, and DyMnOs whereas 
P\\a in Eui_xYxMn03 and Gdo.rTbo.sMnOs are ob- 
served, and then E\\a and E\\c absorptions are expected, 
respectively. However, in reality, the str ong est absorp- 
tions are always observed for E\\a [i,[i,[i3,[ill[l3,[ll,[ii|- 
Moreover, recent experiments in DyMnOs confirm that, 
even when the polarization fiops from P\\c to P\\a by 
the external magnetic fields, the selection rule hardly 
changes [12 1. 

In this paper, we propose a mechanism of absorption 
originated from low-energy spin excitations induced by 
electric fields at terahertz frequencies based on a theory 
proposed for far infrared absorptions in antifcrromag- 
magnon absorptions due to vir- 
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tual electron hoppings are taken into account while we 
neglect the phonon assisted magnon absorption. Wc as- 
sume that antifcrroelectric spin-dependent polarization 
proportional to symmetric spin term Si ■ Sj is domi- 
nant. As observed in Neel ordered state, the conven- 
tional simultaneous two-magnon absorption is possible 
even in cycloidal spin state. However, our results in- 
dicate that novel one-magnon absorption induced by the 
electric component of light is dominant in most of the cy- 
cloidal spin states. Because of the crystal structure, only 
E\\a component produces the one-magnon absorption in 
i^MnOs. We can reproduce some parts of experimen- 
tal observation by considering a one-magnon absorption 
quite well. 

The spin Hamiltonian under the external electric field 
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FIG. 1: (Color online) (a) Cycloidal spin structures, orbital 
ordering pattern on Mn site, and direction of symmetric spin 
dependent polarizations, (b) Non coUinear spin structure on 
a dimer. Applying an electric field along the bond produces 
an effective transverse staggered fields H"^. (c) For E\\a, signs 
of ■ n) are shown on each bond. Spins are oscillated by 
an effective transverse fields shown by (red) arrows, which 
induces one-magnon resonance due to electric fields, (d) For 
E\\b, signs of (£?" ■ II) are shown. The effective fields cancel 
out. 



can be written as H ~ Hq — E ■ P, where Hq is a spin 
Hamiltonian like a Heisenberg Hamiltonian, and P is an 
electric dipole moment which depends on the spin con- 
figurations. In i?Mn03, the electric dipole moment Pij 
associated with a pair of spins on the nearest-neighbor 
bonds have been calculated by a microscopic theory for 
the Mn-O-Mn 180° bond In spin dependent polar- 
izations, the symmetric spin dependent term 



(1) 



is likely dominant 1^, l3| . Since such a kind of polariza- 



tion is realized due to the 



orbital or- 



dering at Mn sites, the directions of 11 are determined as 
shown in Fig.[T](a). For the cycloidal spin structure, the 
local symmetric spin polarization terms Ps are aligned 
antiferroelectrically (see Fig.[l](a)). This additional term 
does not violate the fact that the total fcrroelectricity oc- 
curs solely due to the antisymmetric spin dependent term 
Pa (x Si X Sj. 

We can easily see that one-magnon can be induced 
by the electric component of light for noncoUinear spin 
structures. Let us concentrate on one bond connecting 
site I and 2 (Fig. [1] (b)). When we apply the oscillating 
electric field E'^ along the bond, Hamiltonian of the field 
is -i?" • A = -(^" • n)S^i • ^2. Since Si and ^2 are non- 
coUinear, the spin Si is oscillated by an effective field per- 
pendicular to S'l, H^j_ = (£"^ •n)S'sin6'. The same calcu- 
lation for the spin 5*2 indicates that the effective oscillat- 
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FIG. 2: (Color online) The Heisenberg model for iiMnOa 
with nearest neighbor interactions Ji, antiferromagnetic in- 
teractions for 2nd neighboring b direction J2(> 0), and anti- 
ferromagnetic inter layer interactions Jc(> 0). The direction 
of Ilij are shown by thick (red) arrows. The direction of flij. 
The thin (blue) arrows describe the ground state spin config- 
uration for the be cycloidal state. 



electric component of light can induce a one-magnon res- 
onance caused by effective staggered transverse fields, 
once a noncoUinear ground state is realized. 

The symmetric spin dependent polarization gives us 
the selection rule E\\a straight forwardly. In i?Mn03, the 
propagation vector of cycloidal state is aligned to h direc- 
tions, i.e., spins along a direction are uniform. When we 
apply the electric field E\\a, i?" -11 is uniform (staggered) 
along a (b) direction. As a result, electric fields couple 
with the spin structures. Thus, spins are oscillated with 
the effective staggered fields (Fig. [1] (c)), which induces 
one-magnon resonance. On the other hands, for E\\b, 
E'^ ■ n is staggered (uniform) along a [h) direction. Stag- 
gered -E" • If along a direction docs not couple with the 
uniform spins along a direction, i.e., the effective fields 
cancels out, and, thus, no one-magnon resonance occurs. 

To proceed more quantitative argument in i?Mn03, let 
us consider the frustrated 3-dimcnsional S — 2 Heisen- 
berg Hamiltonian as shown in Fig. [2] as the simplest 
model to reproduce its magnetic behaviors. J2 and Jc are 
fixed to be antiferromagnetic interactions. Although Ji is 
a ferromagnetic interaction in i^MnOa, we also extended 
our calculation to the antiferromagnetic Ji region. In 
i?Mn03, interaction J2 strongly depends R ions. Thus, 
we treat the ratio J2/J1 as a parameters and fixes the 
other parameters to be |Ji| = 8 meV, and Jc/\Ji\ = f.5 
which well reproduces the mag netic behaviors of i?Mn03 
including phase diagrams [2C 



2l| . Due to the frustra- 



ing transverse fields arc staggered: 



-77", . Thus, 



tion, a spiral spin state is a ground state in the param- 
eter range J2/IJ1I > 0.5 for classical spins, and a spiral 
angle 9 is defined as cos 6' = — J1/2J2. An anisotropy 
term {Sf) is added to fix the spiral plane to be he 

{D/\.h\=0.2). 

At low temperatures, the magnetic behaviors of cy- 
cloidal ordered magnets can be described by the lin- 
ear spin wave theory (2^ . After rotations of local axis 
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and Holstein-PrimakofF approximations, the spin Hamil- 
tonian is written by the spin wave annihilation (creation) 
operator (at) [2l| . The spin-dependent polarization is 
also represented by spin wave annihilation and creation 
operators. The polarization ([T]) is written as one- and 
two-magnon processes: 



Ps 



iS^/SN sine n(i)(fc2^)(a^^^ - afc.J 



aka-k-k2 



(2) 



Here n(i) has only a component (see Fig[T](c) and (d)). 
One-magnon absorption can be induced only for the con- 
dition E\\a but two-magnon absorption occurs for both 
E\\a and b. In the one-magnon process, magnons at 
^27r = (27r/a, 0,0) and (0, 27r/6, 0) are induced and each 
of them has the same magnon energy uj2tv Details of 11(2) 
will be reported elsewhere. 

The imaginary part of complex electric polarizability 
tensor at T = 0, which represents an absorption, is ob- 
tained from Kubo formula: ImXaai^) = ^'n^Xalij-^) + 
Imxi^2(w) where 



Imxi^2(w) 



k 

x6{uj - uJk - uj-k-2TT) (4) 



The first term Imxa2(w) represents the one-magnon ab- 
sorption. The position of the absorption peak exists at 
a band edges at k2TT- As a function of qo = 0/{2TT/b), 
one-magnon peak positions are shown in Fig. [31(a). The 
second term corresponds to the absorption due to simul- 
taneous two-magnon absorption process. 

To compare the effects of one- and two-magnon pro- 
cesses for E\\a, integrated intensities defined as 



/ lmxlf!,Kco)du (/3 = 1 or 2) 



(5) 



have been calculated as a function of go, respectively. The 
results are shown in Fig. [3] (b) . Obviously, one-magnon 
absorption is quite strong in a wide range of parameters 
except for near ferromagnetic phase go ~ and antifer- 
romagnetic phase go ~ 1- The effects of two magnon 
absorption are negligible except near the Neel ordered 
state as read off from /f^) //(i) in Fig. [3] (a). The behav- 
iors of the two magnon absorption intensity are consistent 
with the well-known fact that it is observed in antiferro- 
magnets but not in ferromagnets [lil . 17 . [lit . We con- 
clude that for the ferromagnetic interaction Ji (go < 0.5) 
region which iJMnOa belongs to, the two-magnon ab- 
sorption is negligible, contrary to our first expectation 







\ 










• - 






|(2),|(1) . 


t 




0.2 0.4 0.6 0.8 1 

(a) 



FIG. 3: (Color online) (a) The position of the delta function 
of a one-magnon absorption aj27r are shown by (red) squares. 
]Ji| is taken to be 8 meV. The ration of j''^^/ /'^^ are also 
plotted by filled (black) circles, (b) Intensities of one-magnon 
absorption /'^^ shown by (green) circles and two-magnon ab- 
sorption /'^^ shown by (blue) triangles. 



based on the magnon DOS that the absorption originates 
from the conventional two magnon absorption |14| . In a 
view point of the selection rule, dominancy of this one- 
magnon excitation is consistent with the experimental 
results since light absorption occurs only for i?||a and it 
is independent of the direction of the cycloidal spin plane. 

Let us discuss the spectral shapes. As typical exam- 
ples, we compare Imxia'' with imaginary parts of e/i spec- 
tra in DyMnOs and TbMnOa on the assump- 
tion that Ime/.j ~ e" = Imxia"*. We use interactions 
J2I J\ = — 1.2 and J2lJ\ = —0.8 to reproduce cycloidal 
angles in DyMnOs (go ^ 0.36) and TbMn03(go ^ 0.29) 



respectively. Calculated spectra Imxia (arb. unit) to- 



are 



gether with experimental data in Refs. [1^ and 114 
shown in Figs. H] (a) and (b). Delta functions are re- 
placed by the Lorentzian with width 6=1. The disper- 
sions for magnons are shown in Fig. |4] (c). For TbMOs, 
the results are consistent with the experimental data in 
Ref. [2^. As shown in Fig. 3] (b), our results reproduce 
the higher energy peak observed in TbMnOa at a; ~ 8 
mcV. In addition, the shoulder structure in DyMnOs at 
uj ^ h mcV might be interpreted as one-magnon absorp- 
tion as in Fig. |4] (a). However, the origin of the lower 
energy peaks with substantial oscillator strength around 
2 meV is not clear yet. More obviously in DyMnOa, one- 
magnon absorption mechanism alone can not explain the 
whole spectra in i?Mn03. 

Finally, we mention about the effects of anisotropy 
terms. Once the superstructure is induced by an 
anisotropy, one may think that such low energy peaks 
can be generated because a magnon at zone boundary 
mixes with magnons around zone center. However, satel- 
lite peaks induced by superstructures are likely too small 
to reproduce an experimental observation. For example, 
in i?Mn03, a spin fan structure is proposed by includ- 
ing the effects of anisotropy 0, 21 1 , which explains well 



4 








Y.Takahashi m at. ■■ >< « < 











(c)iO| , , , , 1 (d) 

JV1J,I = O.S. q^. = . 

f^Jl\ = O S, 9, = 7 . 

= 7.2, 5, =0 

JVU,I = 1.2. (/,, = / ...^^ 





0.2 0.4 0.6 0.8 1 



FIG. 4: (Color online) (a) Comparison of spectrum shapes 
for DyMnOs between calculation for J2/IJ1I = 1.2 and ex- 
perimental data Im[e/i] extracted from Ref . [l^. (b) Compar- 
ison of spectrum shapes for TbMnOa between calculation for 
J2/I Jil = 0.8 and experimental data Im[e^] extracted from 
Ref. [ij. (c) Dispersion relation of the magnon along 6 di- 
rection as a function of at qa — 0. Circles indicates the 
magnon energy induced by electric field, (c) Spectrum shapes 
for J2/I Ji| = with an anisotropy Dt- 



an elliptical modulation observed experimentally 13, 2^ ■ 
Thus, we examine the effect of easy axis anisotropy along 
b axis -E^b(^')' for J2/IJ1I = -I/V2, (go = 0.25), 
where a fan structure is realized. The results are shown 
in Fig. [4] (d) . The absorption at low energies appears but 
the intensity is too small compared with the observed 
peak in TbMnOs even in the case for a rather strong 
anisotropy, i.e. £)f,/|Ji| = 0.6. In this way, such an ef- 
fect produces a smaller satellites and fails to explain a 
larger peak as in DyMnOs except for the unlikely case 
that modulations of the cycloidal spin structure are ex- 
tremely large. 

So far, it has been believed that the electric compo- 
nent of light can not induce a one-magnon resonance in 
Heisenberg systems through 5^; • Sj, because the sym- 
metric spin polarization ([T]) conserves the total spin mo- 
ment as seen in a two-magnon absorption in Neel ordered 
18| . However, we can now predict that one- 
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magnon resonances can be induced by oscillating electric 
field even in some non-coUinear Heisenberg system. This 
becomes possible because the ground state itself is sym- 
metry broken, which makes a matrix element between 
ground state and one-magnon excitations non-zero. In 
this paper, we have clarified that one-magnon absorp- 
tion is possible for a certain type of cycloidal spin states. 
Moreover, one-magnon resonance might also be observed 
in some other non-colinear spin systems, e.g. canted Neel 
ordering state as we can imagine from Fig. [T] (b). Our 
results indicate the possibility that electromagnons are 
detectable in a wide range of materials. 



In terahertz absorption in i?Mn03, a few mystery are 
still left as open issues. The most challenging and im- 
portant point is to clarify the origin of the lower en- 
ergy absorption around 2meV. Also in absorption around 
5 ~ 8meV, we need a rather large damping factor to re- 
produce experimental results, but the reason of its large- 
ness is not obvious. Such a rather widely spreaded spec- 
tra are also observed even in conventional magnon by the 
inelastic neutron scattering experiments [2^, thus, that 
might be come from the peculiar features of the magnon 
in the cycloidal spin state. The E\\a absorption is ob- 
served even in the high temperature coUinear phase, i.e., 
paraelectric phase 12, 3]. Within the linear spin wave 
theory, we can not treat the coUinear phase properly. 
However, if we assume that the coUinear phase is real- 
ized as a superposition of nearly degenerate ab and be 
cycloidal planes as proposed in Refs.l2d andl2ll. the one- 
magnon absorption due to symmetric term might occur 
even in the coUinear phase within our framework. 
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ment, and N. Kida, M. Mochizuki, Y. Takahashi, T. 
Arima, R. Shimano, and Y. Tokura for fruitful discus- 
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for Scientific Research from the Ministry of Education, 
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Note added — Recently, we became aware of similar 
work by R. Valdcs Aguilar [2^, where they report on 
absorption in TbMnOs and indicate the existence of one- 
magnon absorption considering the effects of orthorhom- 
bic lattice distortions. 
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